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SECTION  1 
INTRODUCTION 


Information  on  the  reflection  of  spherical  blast  waves  is  important  for  explosion 
safety  analysis  and  for  military  applications.  For  example,  one  needs  to  know  the 
peak  pressure  as  a  function  of  distance  from  an  explosion  and  how  it  varies  as  a  func¬ 
tion  of  height  of  burst  (HOB).  Such  information  can  be  conveniently  summarized 
as  isobaric  curves  in  the  height  of  burst-ground  range  plane. 

The  transition  region  between  regular  reflection  and  Mach  reflection  is  especially 
interesting  from  a  gasdynamic  point  of  view.  The  reflection  process  can  lead  to 
significant  pressure  enhancements  in  the  low-pressure  regime  (Ap  <  1  bar).  These 
effects  appear  as  extended  “knees”  in  the  isobaric  HOB  curves.  The  technical 
question  is  the  extent  of  these  “knees”,  and  whether  they  might  be  influenced  by 
read  surface  effects  such  as  boundary  layers  and  surface  roughness. 

Considerable  information  exists  on  the  reflection  of  spherical  blaist  waves  from  pla¬ 
nar  surfaces.  In  the  low-pressure  regime,  HOB  tests  have  been  performed  with 
1-kg  RDX  spheres  (Heilig  et  al.,  1985)  and  with  1-kg  TNT  spheres  (Kingery  and 
Bulmash,  1984;  Petes,  1986).  In  the  high-pressure  regime,  HOB  tests  were  con¬ 
ducted  with  3.6-kg  PBX-9404  spheres  (Carpenter  and  Brode,  1974)  and  with  450- 
kg  PENTOLITE  spheres  (Reisler  et  al.,  1988).  Because  of  charge  size,  these  tests 
were  necessarily  performed  outdoors  and  were,  therefore,  expensive,  small  in  num¬ 
ber,  and  subject  to  weather  variations.  Also,  these  studies  did  not  systematically 
investigate  real-surface  effects  on  the  transition  process. 

This  paper  describes  a  series  of  small-scale  HOB  experiments  developed  at  the 
Emst-Mach  Institut  (Reichenbach  and  Kuhl,  1989).  The  length  scale  of  the 
experiments  is  approximately  1/1000  per  kiloton.  The  advantages  of  such  labora¬ 
tory  experiments  are:  (1)  they  are  inexpensive,  hence  extensive  parametric  studies 
are  feasible  (e.g.,  we  have  conducted  more  than  700  HOB  tests  in  nine  months); 
(2)  weather  variations  are  eliminated  (this  is  especially  important  in  the  low-pressure 
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regime);  (3)  the  roughness  of  the  reflecting  surface  is  easily  controlled;  (4)  high- 
fidelity  photographic  visualization  is  possible. 

The  experiments  axe  described  in  Section  2.  The  results  for  a  hydrodynamically- 
smooth  surface  are  presented  in  Section  3  and  compared  with  other  HE  test  results 
in  Section  4.  The  effects  of  surface  roughness  and  porosity  are  described  in  Section 
5.  The  NP  data  scaled  to  1-KT  conditions  are  shown  in  Section  6,  and  compared 
with  data  and  hydrocode  calculations  in  Section  7.  Our  conclusions  and  recommen¬ 
dations  are  presented  in  Sections  8  and  9. 


2 


SECTION  2 

EXPERIMENT  DESCRIPTION 


The  experiments  were  performed  in  the  cubical  explosion  chamber  (120-cm  by 
120-cm  by  60-cm)  depicted  in  Figure  1.  This  provided  well-controlled  ambient  condi¬ 
tions,  and  protected  personnel  and  equipment  from  the  effects  of  the  explosions.  The 
chamber  was  equipped  with  optical-quality  windows  (69-cm  by 
20-cm  by  2-cm  thick)  for  high-speed  photography.  The  reflecting  plane  was  lo¬ 
cated  on  a  replaceable  roof  section  of  the  chamber.  Three  different  reflecting 
surfaces  were  used:  (1)  a  hydrodynamically-smooth  surface  (Makrolon;  surface 
roughness,  e  <  0.1/zm);  (2)  a  hydrodynamically-rough  surface  (#40  sandpaper: 
e  ~  400  fim  =  5  x  10'3  m/kg1/3  =  1  ft/KT1/3);  (3)  a  porous  rough  surface  (1-cm 
thick  Filtrokelit  ceramic  plate:  e  ~  30  fim  =  3.8  x  10-4  m/kg1/3  =  0.08  ft/KT1/3, 
and  porosity  =  24%  by  volume). 

A  new  explosive  was  developed  for  the  experiments.  It  consisted  of  extremely- 
fine-powdered  NITROPENTA  (NP)  dissolved  in  a  mixture  of  gun  cotton  and  ace¬ 
tone.  To  produce  a  spherical  charge  with  central  ignition,  the  procedure  shown  in 
Figure  2a  was  developed.  Two  electrode  wires  were  centered  in  a  mandrel  holder  (i); 
a  10- /zm  diameter  ignition  wire  was  wrapped  around  them  and  fixed  with  a  silver 
solution.  The  holder  was  dipped  several  times  (ii-iv)  into  the  explosive  mixture, 
which  had  the  consistency  of  honey,  until  an  ellipsoidal  shape  was  achieved  (v). 
After  drying,  the  mandrel  was  placed  in  a  lathe  and  the  charge  was  machined  into  a 
sphere  of  1-cm  diameter.  Then  the  mandrel  holder  was  removed  (vi).  This  method 
allowed  us  to  position  the  ignitor  wire  precisely  at  the  center  of  the  charge.  The 
finished  density  of  the  charge  was  1.40  ±  0.05  g/cm3.  The  finished  charge  mass  was 
approximately  0.5-g;  this  gives  a  length  scale  factor  of  12.6  to  a  kilogram  charge,  or 
854  to  a  kiloton  point  explosion. 

The  surface-burst  charges  were  constructed  by  a  different  method  (Fig.  2b).  A  mold 
was  used  to  cast  0.5-g  hemispherical  charges  of  NITROPENTA  with  a  diameter  of 
1.26  cm.  To  avoid  any  perturbations  on  the  blast  wave,  the  electrode  wires  were 
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connected  to  the  charge  from  below.  To  isolate  the  chamber  from  these  high-voltage 
wires,  a  nonconducting,  disposable  baseplate  of  plexiglass  was  used.  The  electrode 
wires  were  fed  through  two  holes  in  the  baseplate  to  a  small  cavity  provided  for  the 
ignitor  wire.  The  cavity  was  filled  with  a  few  milligrams  of  pure  NITROPENTA 
that  acted  as  a  booster,  and  then  the  charge  was  glued  to  the  baseplate.  Thus, 
the  ignitor  wire  was  positioned  at  the  center  of  the  base  of  the  hemisphere.  The 
baseplate  was  mounted  in  a  disposable  Delrin  cylinder  that  could  deform  during 
the  explosion.  The  entire  assembly  was  emplaced  in  the  chamber  by  means  of  a 
mounting  bracket  (Fig.  2c)  that  shock-isolated  the  reflecting  surface  from  the  stress 
waves  produced  by  the  charge.  The  disadvantage  of  this  system  was  that  two  holes 
in  the  plexiglass  baseplate  created  a  preferential  relief  path  for  the  stresses  under  the 
charge.  This  permitted  significant  cratering  of  the  plexiglass,  with  a  concomitant 
loss  of  approximately  25  percent  of  the  energy  from  the  blast  wave  (see  Section  3). 

The  electrode  wires  served  as  the  charge  support  system  for  the  height-of-burst 
tests.  The  charge  position  was  controlled  within  ±0.1  mm  by  machined  gauge 
blocks.  The  charge  was  fired  by  dumping  about  50  joules  energy  from  a  10  kV 
capacitor  bank  into  the  fine  ignitor  wire.  This  created  an  over-driven  detonation 
wave  that  caused  a  complete  and  repeatable  detonation  of  the  entire  charge  mass. 
The  timing  and  firing  was  controlled  by  a  computer  that  also  synchronized  the 
pressure  recording  and  photographic  systems. 

The  single-frame  photography  system  consisted  of  a  single-spark  light  source 
(300-ns  duration),  two  spherical  mirrors  (64-cm  diameter,  400-cm  focal  length) 
which  produced  a  parallel  light  beam  through  the  test  section,  and  an  open-shutter 
camera  (Fig.  3).  This  system  gave  high-quality  shadow  pictures.  By  placing  a  ra¬ 
zor  blade  edge  at  the  second  focus,  excellent  schlieren  photographs  could  also  be 
produced. 

Static  pressure  histories  were  measured  at  13  ground  ranges  along  the  reflecting 
surface.  Kistler  transducers  (type  603B)  with  charge  amplifiers  (type  5001)  were 
used.  The  sensitive  area  of  the  gauges  was  about  3  mm  in  diameter,  thus,  wave¬ 
form  features  that  were  smaller  than  3  m/KT1/3  or  8  ms/KT1/3  could  not  be 
resolved  accurately.  The  data  were  recorded  on  a  LeCroy  digital  transient  recorder 
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(type  8837F)  which  had  a  sampling  frequency  of  10  MHz.  The  gauge  output  was 
filtered  at  280  kHz;  under  normal-shock  reflection  conditions,  this  system  had  a 
rise  time  of  1  ps  without  overshoot  in  the  record.  Special  precautions  were  taken 
to  shock-isolate  each  pressures  gauge  from  the  blast-induced  vibrations  in  the  re¬ 
flecting  plane,  e.g.,  Figure  4  shows  that  no  metal  contact  was  allowed  between  the 
pressure  gauge  support  and  the  roof. 

The  pressure  measurements  were  systematically  analyzed  by  the  following  proce¬ 
dure.  For  each  experiment,  the  measured  peak-pressures  (i.e.,  without  extrapola¬ 
tion)  were  tabulated  as  a  function  of  ground  range  and  scaled  to  sea- level  conditions. 
Then  the  data  were  fit  with  a  least-squares  computer  program  to  produce  an  ana¬ 
lytic  fit  to  the  peak  pressure- range  curve  for  each  height  of  burst.  These  analytic 
functions  were  fed  into  the  HOB  program  that  calculated  the  iso-overpressure  ranges 
and  plotted  the  data  as  height-of-burst  curves. 

Experiments  were  performed  for  22  different  heights  of  burst,  starting  with 
0.5-g  NP  hemisphere  surface  bursts,  and  extending  to  heights  of  burst  of  500  mm 
(6.3  m/kg1/3  ~  1400  ft/kT1/3)  relevant  to  the  low-pressure  regime.  First,  a  baseline 
set  of  experiments  was  conducted  using  the  Makrolon  plate  as  a  reflecting  plane. 
This  data  served  as  a  benchmark  control  case  for  a  hydrodynamically-smooth  sur¬ 
face.  Then  the  test  series  was  repeated  over  a  rough  surface  (#40  sandpaper)  and 
over  a  porous  Filtrokelit  plate — to  investigate  how  boundary  layer  momentum  losses 
and  mass  losses  affect  the  HOB  curves. 
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SECTION  3 

SMOOTH-SURFACE  RESULTS 


Figure  5  presents  shadow  photographs  of  the  blast  wave  reflection  from  the  Makrolon 
surface  in  the  low  pressure  regime  (HOB  =  303  mm  =  3.8  m/kg1/3).  Figure  5a 
depicts  the  regular  reflection  of  the  incident  shock  I.  Also  shown  is  the  backward- 
facing  shock  Ij  that  is  characteristic  of  HE-driven  blast  waves,  and  its  reflection 
from  the  surface  (shock  Ri).  Figure  5b  depicts  the  shock  structure  after  transition 
to  Mach  reflection;  the  Mach  stem  NR  and  the  slipline  SL  are  clearly  visible.  As 
will  be  discussed  in  Section  8,  this  is  actually  a  Neumann  reflection  shock  structure 
and  not  a  classical  3-shock  Mach  configuration.  Shock  R'  denotes  the  reflection  of 
shock  R  from  the  chamber  floor  (recall  that  the  tests  were  performed  upside  down). 
Secondary  images  of  the  shock  structures  are  caused  by  the  curved-shock  effects. 

Figure  6  presents  the  peak  overpressures  (Apr)  measured  at  ground  zero  as  a  func¬ 
tion  of  height  of  burst.  Using  normal  shock  reflection  factors,  this  curve  was  “un¬ 
reflected”  to  create  the  incident  pressure-range  curve  ( Api )  for  the  NITROPENTA 
blast  wave  (see  Eq.  (10)  for  an  analytic  fit  to  this  free-air  curve). 

Figure  7  depicts  the  peak-pressure  height-of-burst  curves  for  the  hydrodynamically- 
smooth  surface  in  the  low-pressure  regime  (A p  =  1  to  15  psi).  The  circles  denote 
the  measured  points,  while  the  symbols  “+”  and  ”  represent  the  -f  1  a  and  — l<x 
uncertainty  bands  on  each  pressure-range  curve  (typically  about  ±  1  percent  in 
range).  The  data  scatter  is  extremely  small,  and  there  are  so  many  data  points 
that  it  is  very  easy  to  draw  the  HOB  curves  through  the  data.  The  long  dashed  line 
denotes  the  limit  of  measurement  region;  the  2,  1.5  and  1  psi  curves  were  obtained 
by  extrapolation;  a  larger  chamber  is  needed  to  measure  these  low  pressures.  The 
short  dashed  dotted  lines  denote  the  HOB  curves  in  the  regular  reflection  region — as 
calculated  from  the  von  Neumann  two-shock  theory  and  the  incident  free-air  curve 
of  Figure  6.  The  experimental  HOB  curves  compare  well  with  the  two-shock  theory 
for  about  the  first  half  of  the  regular  reflection  region,  but  in  the  transition  region 
between  regular  and  Mach  reflection,  the  experimental  curves  fail  systematically 
below  the  theoretical  curves  for  pressures  less  than  15  psi. 
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For  the  surface  burst  case  (HOB  =  0),  there  axe  always  two  data  points  for  each 
pressure  level:  the  circle  (the  left-most  point)  denotes  the  measured  pressure  from 
the  0.5-g  hemisphere  experiments;  the  star  (the  right-most  point)  corresponds  to  an 
ideal  surface  burst.  It  was  calculated  from  the  NP  free-air  curve  by  multiplying  the 
range  by  21/3  (i.e.,  an  ideal-surface  approximation).  The  measured  surface  burst 
pressures  occur  at  a  ground  range  of  about  90  percent  of  the  ideal-surface  case,  i.e., 
GRsb/GRisb  =  0.90;  cubing  this  gives  the  explosion  energy  ratio:  Esb/Eisb  — 
( GRsb/GRisb )3  =  0.751.  Thus  only  about  75  percent  of  the  chemical  energy 
for  the  hemispheric  surface-burst  charge  was  available  to  drive  the  blast  wave,  and 
25  percent  of  the  energy  was  lost  to  the  cratering  process. 

Figure  8  presents  the  height-of-burst  curves  for  the  hydrodynamically-smooth  sur¬ 
face  in  the  intermediate-pressure  regime  (A p  =  8  to  200  psi).  The  notation  is  the 
same  as  in  the  previous  figure.  Again  the  data  scatter  is  extremely  small,  and  the 
HOB  curves  are  very  smooth.  Again  we  see  that  for  low  pressures,  the  experimen¬ 
tal  HOB  curves  Eire  below  the  theoretical  curves  in  the  transition  region,  while  for 
higher  pressures  the  data  is  in  agreement  with  the  theory  throughout  most  of  the 
regular  reflection  regime. 

The  charge  center  for  the  lowest  four  HOB  experiments  was  located  at  heights  of 
0.5,  1,  1.5  and  2  cm;  these  correspond  to  a  surface  tangent  case  and  to  heights  of 
burst  of  1,  1.5  and  2  charge  diameters,  respectively.  Because  these  charges  were  so 
close  to  the  reflecting  plane,  a  shaped-charge  or  piston  effect  was  created  by  the 
dense  detonation  products  jetting  along  the  surface.  This  piston  effect  caused  the 
shock  to  bulge  outward  near  the  surface,  and  distorted  the  pressure-range  curve. 
This  “HE  squish  effect”  appears  as  a  bulge  in  the  HOB  curves  for  low  burst  heights 
HOB  <  3.3  cm  =  0.4  m/kg1/3).  Figure  8  shows  that  the  squish  effect  is  most 
prominent  at  high  pressures,  and  decreases  at  larger  ground  ranges. 

Figure  9  presents  the  height-of-burst  curves  for  the  hydrodynamically-smooth  sur¬ 
face  in  high-pressure  regime  (Ap  =  100  to  1000  psi).  The  notation  is  the  same 
as  in  the  previous  figures.  The  “HE  squish  effect”  is  especially  prominent  in  this 
figure.  At  these  high  pressures,  the  data  are  consistent  with  the  two-shock  theory 
in  the  regular  reflection  regime.  The  200  psi  data  point  at  HOB  =  1.2  m/kg1/3  is 


7 


especially  supportive  of  this  statement.  However,  more  experiments  at  intermediate 
heights  of  burst  are  needed  to  accurately  define  the  knees  in  the  NP  height -of-burst 
curves  in  the  high-pressure  regime. 
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SECTION  4 

COMPARISONS  WITH  OTHER  HE  DATA 


In  this  section  we  shall  compare  the  NP  height-of-burst  curves  with  data  taken  on  a 
variety  of  HE  field  tests  performed  at  much  larger  scale.  We  used  the  same  method 
of  analysis  for  the  HE  field  tests  that  we  had  used  for  the  NP  experiments;  namely, 
the  data  were  cube-root  scaled  to  sea-level  conditions  and  then  peak-  pressure  versus 
ground  range  data  were  curve  fit  using  a  least-  squares  program.  Then  the  analytic 
fits  were  interrogated  to  find  the  iso-overpressure  ground  ranges  and  compared  with 
the  NP  height-of-burst  curves. 

The  first  set  of  HE  data  that  we  analyzed  is  listed  in  Table  1.  The  height-of- 
burst  data  came  from  the  DIPOLE  WEST  tests  using  450-kg  charges  of  TNT  and 
the  MIGHTY  MACH  tests  employing  primarily  29-kg  charges  of  PENTOLITE; 
both  were  performed  over  concrete  pads.  The  surface-burst  data  came  from  tests 
employing  much  larger  stacked  TNT  charges.  In  those  cases,  a  considerable  fraction 
of  the  energy  was  lost  in  the  cratering  process.  Because  the  scatter  was  typically 
quite  large  in  this  data  set,  we  did  not  take  into  account  the  energy  differences 
between  TNT  and  NITROPENTA. 

Figures  10  and  11  compare  this  HE  data  set  with  our  smooth-  surface  NP  curves. 
At  low  pressures  (Fig.  10),  the  HE  data  is  in  excellent  agreement  with  our  NP 
curves — not  only  in  Mach  reflection  region  (i.e.,  below  the  knees)  but  also  in  the 
transition  region  and  at  the  knees  of  the  curves.  Most  notable  are  the  pressure 
levels  of  Apr  =  1,  1.5,  2  and  4  psi.  In  the  intermediate  pressure  regime  (Fig.  11), 
the  HE  data  are  also  consistent  with  the  NP  curves  but  the  HE  data  have  much 
more  scatter,  especially  for  pressures  of  Apr  =  8  to  15  psi.  At  somewhat  higher 
pressures  (  Apr  =  30  to  70  psi),  the  HE  data  lie  just  inside  the  NP  curves  but  follow 
their  general  shape.  In  fact,  the  NP  curves  are  closer  to  the  two-shock  theory, 
especially  near  transition,  than  the  large-charge  data — perhaps  because  we  used 
a  hydrodynamically-smooth  surface.  Thus,  these  figures  show  that  increasing  the 
explosion  length  scale  by  a  factor  of  100  does  not  raise  the  knees  of  the  height-of- 
burst  curves  in  the  low-pressure  transition  region  (because  of  the  overall  system 


9 


Table  1.  HE  test  data  sources. 


w 

(lb) 

HOB 

(m/kg1/3) 

Charge 

Event 

0 

TNT 

FLATTOP  II,  FLATTOP  III,  DPL  3,  MIDDLEGUST  I1  1 

.05 

TNT 

DIAL  PACK,  DPL  6A,  DPL  6,  PRAIRIE  FLAT,  MIXED  CO.1 

.10 

TNT 

MIDDLEGUST  II,  MINEUNDER1 

~  1000 

.205 

TNT 

Shot  10  (1969  SERIES)1 

-  1000 

.28 

TNT 

Shot  9  (1969  SERIES)1 

-  40, 000 

.33 

TNT 

Shot  1A  (DISTANT  PLAIN)1 

-  1000 

.58 

TNT 

Shot  8  (1969  SERIES)1 

~  1000 

.91 

TNT 

Shot  1,  Shot  2  (1969  SERIES)1 

-40,000 

.97 

TNT 

Shot  1  (DISTANT  PLAIN)1 

-  1000 

1.102 

TNT 

Shot  4  (1969  SERIES)1 

~  1000 

1.38 

TNT 

Shot  3  (1969  SERIES)1 

1040 

1.73 

TNT 

Shot  17  (DIPOLE  WEST)1'2 

-  1000 

1.77 

TNT 

Shot  5  (1969  SERIES)1 

-  1000 

2.26 

TNT 

Shot  6  (1969  SERIES)1 

-  1040 

2.30 

TNT 

Shot  18  (DIPOLE  WEST)12 

-  1000 

2.70 

TNT 

Shot  7  (1969  SERIES)1 

1040 

3.42 

TNT 

Shot  24  (DIPOLE  WEST)12 

1040 

3.45 

TNT 

Shot  19  (DIPOLE  WEST)1-2 

1040 

4.60 

TNT 

Shot  20  (DIPOLE  WEST)1-2 

1040 

4.60 

TNT 

Shot  23  (DIPOLE  WEST)12 

1086 

5.03 

TNT 

Shot  6  (DIPOLE  WEST)1 

1040 

5.49 

TNT 

Shot  21  (DIPOLE  WEST)1'2 

1040 

5.52 

TNT 

Shot  22  (DIPOLE  WEST)1-2 

63 

5.7 

PENTOLITE 

Shot  1  (MIGHTY  MACH  II)2 

218 

7.7 

PENTOLITE 

Shot  1  (MIGHTY  MACH  IV)2 

64 

7.72 

PENTOLITE 

Shot  1  (MIGHTY  MACH  III)2 

64 

7.73 

PENTOLITE 

Shot  2  (MIGHTY  MACH  II)2 

64 

7.77 

PENTOLITE 

Shot  1A  (MIGHTY  MACH  III)2 

217 

9.62 

PENTOLITE 

Shot  2  (MIGHTY  MACH  IV)2 

64 

9.65 

PENTOLITE 

Shot  2  (MIGHTY  MACH  III)2 

Superscript:  1  denotes  Carpenter  (1978)  dataset 


2  denotes  Sauer  (1990)  data  set 
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response).  We  infer  from  this  comparison  that  our  peak-pressure  measurements  are 
as  accurate  as  those  made  on  the  large  HE  tests. 

These  observations  are  supported  by  more  accurate  HE  data  presented  in 
Figure  12.  Shown  there  are  results  of  height-of-burst  tests  performed  at  the  Ernst- 
Mach- Institute  test  site  at  Wintersweiler  with  RDX  spheres  (Heilig  et  al.,  1985; 
Scheklinski-Gluck,  1990).  Charge  masses  of  1-kg,  0.128- kg  and  0.016- kg  were  used. 
We  note  that  the  RDX  tests  used  the  same  pressure  gauges  (KISTLER  603B)  and 
recording  system  that  were  used  in  the  NP  experiments.  Although  not  quite  as 
smooth  and  consistent,  the  RDX  data  are  in  agreement  with  the  NP  curves.  Most 
notably,  the  RDX  data  also  suggests  that  the  height-of-burst  curves  begin  to  devi¬ 
ate  from  the  two-shock  theory  well  before  the  detachment  angle  (see,  for  example, 
the  curves  for  Apr  =  8,  10  and  15  psi).  At  intermediate  pressures  (e.g.,  Apr  =  15, 
20  and  30  psi),  the  RDX  curves  have  slightly  more  pronounced  knees.  This  is  a 
subtle  but  discernible  trend  that  may  have  been  caused  by  enhanced  resolution  of 
the  peak-pressure  measurements  on  the  RDX  experiments.  Thus,  increasing  the 
explosion  length  scale  by  a  factor  of  10,  and  using  the  same  pressure  measurement 
system  (thereby  increasing  the  temporal  resolution  by  about  10),  has  little  effect  on 
the  shape  of  the  height-of-burst  curves. 

To  check  the  NP  curves  in  the  high-pressure  regime  (100  to  1000  psi)  we  made 
comparisons  with  data  taken  on  tests  with  3.6-kg  spheres  of  PBX-9404  (Carpenter 
and  Brode,  1974).  Figure  13  shows  that  the  NP  data  are  in  excellent  agreement 
with  the  PBX  curves — especially  for  Apr  =  100,  150  and  200  psi — both  in  the 
knee  region  near  transition  and  in  the  HE  squish  region  near  the  surface.  For 
pressures  greater  than  about  150  psi,  the  shock  structure  transitions  to  a  double- 
Mach  configuration  when  the  detachment  angle  is  reached,  and  double-  peaked 
waveforms  are  created.  Smaller  pressure  gauges  are  required  to  measure  the  details 
of  such  waveforms  accurately. 

In  summary,  we  find  that  blast  wave  generated  by  the  0.5-kg  NP  charges  is  consistent 
and  reproducible,  and  provides  an  excellent  laboratory  simulation  of  the  HE-driven 
blast  waves  employed  in  large-scale  field  tests.  The  accuracy  of  the  pressure  mea¬ 
surements  made  in  our  NP  experiments  was  similar  to  the  accuracy  of  pressure 
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measurements  made  on  HE  field  tests  (e.gM  the  HE  data  always  lies  inside  or  on 
the  NP  height-of- burst  curves  for  a  given  pressure  level).  Comparisons  with  the 
RDX  data  set  indicates  that  knees  of  the  curves  could  be  somewhat  more  pointed. 
Experiments  with  smaller-diameter  pressure  gauges  should  be  made  to  check  this 
trend,  and  to  more  fully  define  the  NP  height-of-  burst  curves  in  the  high-pressure 
regime. 
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SECTION  5 

BOUNDARY  LAYER  EFFECTS 


In  this  section  we  explore  how  boundary  layer  effects  change  the  height-of-burst 
curves.  Rough-wall  boundary  layer  effects  were  evaluated  by  comparing  the  results 
of  NP  experiments  conducted  over  #40  sandpaper  surface  (labeled  “R”)  with  the 
smooth-surface  curves,  as  shown  in  Figure  14.  The  large  surface  roughness  (e  = 
400 nm  =  0.005  m/kg1/3  =  1  ft/KT1/3)  used  here  caused  no  significant  effect  on 
the  peak  pressures  in  the  regular  reflection  region  or  in  the  Mach  reflection  region 
below  the  knees  of  the  curves.  Surface  roughness  clipped  or  rounded  off  the  sharp 
peaks  in  the  transition  region;  this  caused  a  small  but  systematic  depression  of  the 
curves  in  the  knee  region. 

The  effects  of  mass  loss  due  to  surface  porosity  were  evaluated  by  comparing  the 
results  of  the  NP  experiments  performed  over  the  Filtrokelit  surface  (labeled  “P”) 
with  the  smooth- surface  curves,  as  also  shown  in  Figure  14.*  Porosity  caused  more 
dramatic  effects  than  surface  roughness,  especially  in  the  knee  region  of  the  curves, 
and  this  effect  increased  at  the  lower  pressures. 

The  above  comparisons  demonstrate  that  the  surface  boundary  condition  has  a 
measurable  effect  on  the  the  shape  of  knees  of  the  height-of-burst  curves.  The 
smooth- surface  curves  are  the  highest,  followed  by  the  rough-surface  curves,  and 
the  porous- surface  curves  are  the  lowest.  By  induction,  then,  one  would  expect  that 
the  inviscid-surface  curves  would  be  even  higher  than  the  smooth-surface  curves. 
Such  effects  are  consistent  with  boundary  layer  theory.  According  to  the  momentum 
integral  for  a  boundary  layer  behind  a  shock,  momentum  loss  due  to  wall  drag  and 
mass  loss  due  to  porosity  both  have  the  same  negative  boundary  layer  displacement 
effect  (i.e.,  a  weak  rarefaction  effect)  on  the  near-surface  flow.  Apparently  the 
porosity  had  a  larger  effect  than  the  surface  roughness.  Height-of-burst  experiments 
over  rough  surfaces  should  be  extended  to  lower  pressures. 


*Comparison  of  height  of  burst  curves  for  other  blast  wave  parameters  (e.g.,  impulse, 
shock  arrival  time  and  positive  phase  duration)  can  be  found  in  the  Appendix. 
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From  photography  on  the  DIPOLE  WEST  tests  we  know  that  the  triple  point 
trajectory  is  suppressed  on  a  concrete  surface  in  comparison  to  an  ideal  reflecting 
plane  (formed  in  those  tests  by  the  symmetry  plane  between  two  simultaneous 
bursts).  This  suppression  is  no  doubt  caused  by  the  boundary  layer  along  the 
surface.  Note  that  even  our  hydrodynamically-smooth  surface  created  a  boundary 
layer.  We  infer  from  this  that  the  triple  point  trajectory  is  suppressed  even  on  such 
a  smooth  surface  (relative  to  the  triple  point  trajectory  over  an  inviscidly-ideal 
surface);  surface  roughness  or  porosity  will  further  enhance  this  effect.  Apparently, 
suppression  of  the  triple  point  trajectory  erodes  the  peak  pressures  and  pulls  in  the 
knees  of  the  height-of- burst  curves,  as  demonstrated  in  Figure  14.  Hence,  boundary 
layer  effects  seem  to  be  responsible  for  transition  well  before  the  detachment  angle 
is  achieved.  But  because  the  boundary  layer  suppresses  the  triple  point  trajectory, 
the  Mach  stem  is  not  detectable  in  photography  until  well  after  the  detachment 
angle  is  reached.  This  will  be  discussed  more  fully  in  Section  7. 

Further  evidence  of  boundary  layer  effects  is  provided  by  the  fact  that  the  height -of- 
burst  curves,  when  extrapolated  to  the  surface,  do  not  agree  with  the  ideal  surface 
burst  points — even  when  a  hydrodynamically-smooth  surface  is  used  as  a  reflecting 
plane.  This  effect  is  demonstrated  in  Figure  15,  which  depicts  the  peak  pressure 
versus  range  curves  for  the  ideal-surface  case  (label  I  SB)  and  the  extrapolated 
surface  burst  case  (label  ESB).  Clearly,  the  extrapolated  surface  burst  case  decays 
more  rapidly  than  the  ideal-surface  case.  This  can  also  be  seen  in  the  exponent  of 
the  curve  fits  to  the  data; 


&Pisb  =  9.59/GH1734  (1) 

A  pesb  =  10.96/GR1850  (2) 

where  [Ap]  =  bars,  [GR]  =  m/kg1/3  and  3.5  <  GR  (m/kg1//3)  <  10.  For  convenience, 
the  above  relations  may  be  written  in  a  logarithmic  form: 

log  A pjsb  =  a1-bi  log (GRISb)  (3) 

log  Apesb  =  a2-62  log {GResb)  (4) 

where  a 1  =  0.982,  6j  =  1.734,  02  =  1.04  and  62  =  1.85. 
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Let  us  define  the  percent  decrement  in  ground  range  to  the  same  overpressure: 

6  =  100  (GRisb  —  GResb)/GRisb  (5) 


Using  the  pressure-range  curves  (Eqs.  3  and  4),  one  can  show  that 

6  =  230(ai  -  a2  )/b2  +  230(  1  —  &1/62)  log  GRISB  (6) 

or,  after  inserting  the  constants,  one  finds 

6  =  7.2  +  14.6  log(GR/sB )  (3  <  GR  <  10)  (7) 


In  other  words,  the  ground-range  decrement  increases  logarithmically  with  range 
(see  Fig.  16).  It  is  plausible  that  this  increase  with  range  is  caused  by  boundary 
layer  momentum  and  energy  losses  because  the  surface  area  in  contact  with  the 
hemispherical  blast  wave  increases  as  the  square  of  the  shock  radius. 

One  can  try  to  remove  such  boundary  layer  effects  from  the  NITROPENTA  data  by 
increasing  the  ground  range  of  each  data  point  for  the  smooth-surface  case,  accord¬ 
ing  to  Eq.  7.  This  then  defines  our  approximation  to  the  ideal-surface  height -of-burst 
curves  (symbol  I)  shown  in  Figures  17  and  18.  The  boundary  layer  corrections  be¬ 
gin  at  a  pressure  of  A p  =  15  psi  and  the  corrections  increase  logarithmically  with 
range.  They  extend  the  knees  somewhat,  but  the  ideal-surface  curves  still  remain 
below  the  two-shock  theory  curves  near  transition.  Clearly  this  correction  cannot 
eliminate  the  triple  point  suppression,  and  therefore,  it  cannot  fully  remove  the 
boundary  layer  effects  on  transition. 

In  summary,  we  believe  that  boundary  layer  effects  cause  the  knees  in  the  HOB 
curves  to  be  pulled  in  and  down  in  the  low-pressure  regime.  They  probably  trig¬ 
ger  an  early  transition,  and  suppress  the  triple  point  trajectory.  Apparently  this 
happens  even  for  a  hydrodynamically-smooth  surface.  Comparisons  with  inviscid 
hydrocode  simulations  of  the  reflection  of  the  NP  blast  wave  could  help  verify  this 
assertion. 

An  additional  boundary  effect  that  is  worth  noting.  Height-of-burst  tests  have  been 
conducted  over  snow  and  concrete  surfaces  (Wisotski,  1990).  Typically,  1-lb  TNT 
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charges  were  detonated  over  a  deep  snow  layer  (18  to  127  cm).  Pressure  waveforms 
were  measured  at  the  snow  surface.  We  scaled  the  peak -pressure  measurements 
to  sea-level  conditions  and  plotted  the  peak-pressure  versus  range  curves.  These 
were  used  to  construct  height-of-burst  curves,  which  are  compared  with  our  NP 
curves  in  Figures  19  and  20.  These  figures  show  that  the  HOB  curves  for  snow  are 
dramatically  contracted  relative  to  the  NP  curves  (or  the  concrete-surface  curves) 
for  peak  pressures  of  30  psi  to  3  psi.  The  blast  wave  compresses  the  snow  surface 
which  changes  the  effective  angle  between  the  incident  shock  and  the  reflecting 
plane,  and  thus  creates  a  local  rarefaction  wave  near  the  surface.  Also,  the  blast 
wave  performs  work  in  compressing  the  snow,  and  thus  the  blast  wave  loses  energy 
to  the  snow.  This  dramatically  changes  the  height-of-burst  curves.  No  doubt  this 
effect  depends  on  the  depth  of  the  snow  layer.  For  more  realistically  scaled  depths 
of  the  snow  layer,  this  effect  may  be  considerably  reduced. 
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SECTION  6 

SCALING  TO  POINT  EXPLOSIONS 


In  this  section  we  will  scale  the  NITROPENTA  height-of-burst  curves  to  conditions 
appropriate  to  large-scale  point  explosions,  and  compare  these  results  with  other 
data  and  hydrocode  simulations. 


The  fundamental  input  to  the  scaling  procedure  is  the  free-air  pressure  versus  radius 
curve  for  a  1  kiloton*  point  explosion  (Smiley  et  al.,  1981)  depicted  in  Figure  21. 
The  analytic  fit  to  this  curve  is  as  follows: 

{  5.885  -  2.6  log  Rn  for  67  <  Rn  <  141 

7.147  -  4.279  log  Rn  +  0.51(logRn)2  for  141  <  Rn  <  1106  (8) 

2.613  -  1.241  log  Rn  for  1106  <  R„  <  2926 

where  [Ap]  =  bars  and  [R]  =  m/KT1/3.  Next,  we  calculate  a  local  scale  function, 
SF,  which  is  defined  as  the  point-source  radius  to  a  given  pressure  divided  by  the 
NITROPENTA  radius  to  the  same  pressure: 


SF  =  Rn/Rup 


(9) 


Here,  Rn  comes  from  Figure  21  (or  Eq.  8)  and  Rnp  comes  from  the  NITROPENTA 
free-air  curve  (Fig.  6),  or  equivalently: 


(  1.047  —  2.289  log  Rnp 


for  0.9  <  Rnp  <  2.45 


log  A pNP  = 


{  0.954  -  2.253 log  Rnp  +  0.427(log  Rnp)2 


for  2.45  <  Rnp  <  9.8 


l  0.462  -  1.334  log  Rnp 


for  9.8  <  Rnp  <  28 

(10) 


*Recall  that  1-KT  =  1012  calories. 
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where  [Ap^rp]  =  bars  and  [R/vp]  =  m/kg1/3.  As  shown  in  Figure  22,  this  scale 
function  varies  with  distance: 


SF  =  f(RNP )  = 


73.18  -  28.675log(RNp)  for  RNP  <  1.54 
67.75  for  RNp>  1-54 


(11) 


The  scale  function  starts  with  a  value  of  74.4  at  high  pressures,  decreases  to  a  value 
of  67.75  at  about  60  psi,  and  then  remains  constant  (total  variation  of  about  10%). 
To  provide  a  consistency  check,  this  scale  function  was  applied  to  the  data  points 
for  the  NP  free- air  curve  (Fig.  6),  and  compared  with  the  point-source  free-air  curve 
in  Figure  21.  The  agreement  is  excellent  (deviations  of  less  than  one  percent),  so 
clearly  the  scale  function  (Eq.  11)  accurately  transforms  NP  ranges  to  point-source 
ranges. 


With  the  incident  blast  wave  properties  properly  scaled,  an  inviscid  (i.e.,  cube- 
root  of  charge  energy)  scaling  procedure  can  now  be  defined  for  the  reflected  pres¬ 
sure.  One  assumes  that  the  reflected  pressure  values  depend  on  only  two  param¬ 
eters:  the  incident  shock  pressure  (which  is  only  a  function  of  the  scaled  shock 
radius)  and  the  incident  shock  angle  a  with  respect  to  the  reflecting  surface  (a  = 
tan-1  {GR/ HOB}).  Thus  the  scaling  procedure  reduces  to  the  following  relations: 

R„  =  Rnp  ■  SF(Rnp)  (12) 

<*n  =  ONP  (13) 

This  inviscid  scaling  preserves  the  geometric  similarity  between  the  two  blast  waves. 
Hence,  both  the  ground  range  and  height  of  burst  of  any  data  point  i  on  the 
NITROPENTA  HOB  curves  may  be  scaled  to  an  equivalent  point  for  a  point-source 
blast  wave,  according  to: 


GRn(i)  =  GRNP(i )  •  5F[i2/VP(0]  (14) 

HOBn(i )  =  HOBNP{i)  ■  SF[RNP(i)\  (15) 

where 

Rnp( i)  =  \/GRNp{i )2  +  HOBnp(i)2 
Clearly,  this  scaling  neglects  real-surface  effects  such  as  boundary  layers. 
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Figures  23  and  24  present  the  height-of-burst  curves  for  1-KT  point  explosions, 
constructed  from  Figures  17  and  18  by  means  of  the  aforementioned  inviscid  scal¬ 
ing  procedure.  There  are  two  curves  for  each  pressure  level:  the  curve  with  the 
“O"  symbols  represent  the  smooth-surface  case,  and  the  curve  with  the  sym¬ 
bols  represent  our  approximation  to  the  ideal-surface  case  (i.e.,  where  the  boundary 
layer  decrement  effects  have  been  removed  before  scaling).  Even  after  applying  the 
scale  function,  the  curves  remain  smooth  and  self-consistent. 

In  closing,  it  is  worth  pointing  out  that  the  cube-root  scaling  procedure  employed 
here  is  based  on  inviscid  similarity  theory  for  blast  waves  (Sedov,  1959).  It  is 
theoretically  valid  only  if  the  two  blast  waves  can  be  matched  with  a  single  scale 
constant,  for  example,  scaling  results  of  an  x-kg  charge  to  those  of  a  y-kg  charge 
of  the  same  explosive ,  or  scaling  results  of  an  x-KT  point  explosion  to  those  of  a 
y-KT  point  explosion.  Thus,  the  variable  scale  function  SF  used  here  has  no  firm 
theoretical  basis.  In  fact,  similarity  scaling  does  not  even  work  perfectly  for  all  HE 
blast  wave  parameters;  for  example,  the  cube-root-scaled  positive  phase  duration  of 
TNT  blast  waves  is  different  from  those  of  more  energetic  explosives  such  as  PBX- 
9404  (because  TNT  afterburns  significantly  and  PBX  does  not).  Fortunately,  the 
pressure-range  curve  for  most  explosives  is  somewhat  insensitive  to  these  effects  at 
low  pressures  and  their  pressure- range  curves  can  be  collapsed  by  scaling  the  data 
with  the  cube-root  of  the  yield. 

Besides  this  pedantic  point,  there  are  a  number  of  physical  effects  that  cannot  be 
scaled.  For  example,  when  the  reflected  shock  passes  through  the  fireball,  it  heats 
the  gas  due  to  shock  compression  and  changes  its  entropy.  The  irreversible  part  of 
this  local  energy  change  is  then  not  available  to  drive  the  blast  wave.  Now.  an  HE 
fireball  is  10  to  100  times  more  dense  than  the  fireball  of  a  point  explosion,  so  a 
different  fraction  of  the  blast  wave  mass  (and  thus  energy)  will  be  affected  in  the  two 
examples.  Also,  point  explosions  are  not  constant  energy  blast  waves  (because  of 
radiation  losses  from  the  high-temperature  fireball),  whereas  HE  explosions  produce 
blast  waves  of  constant  or  increasing  energy  (due  to  afterburning).  In  addition,  the 
density  gradient  of  the  atmosphere  can  refract  the  blast  wave  of  large-yield  point 
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explosions,  and  hydrocode  calculations  have  shown  that  this  can  affect  the  height- 
of-burst  curves  for  Apn  <  2  psi.  Finally,  boundary  layer  effects — especially  in  the 
transition  region — cannot  be  scaled. 

In  summary,  there  is  no  scaling  theory  that  properly  accounts  for  all  the  physical 
effects.  Hence,  the  variable  scale  function  SF  employed  here  represents  the  best 
available  engineering  approximation  for  comparing  HE  results  with  those  of  point 
explosions. 
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SECTION  7 

COMPARISON  WITH  POINT  EXPLOSION 
DATA  AND  CALCULATIONS 


Next  we  shall  compare  our  NP  curves  with  the  results  of  inviscid  hydrocode  sim¬ 
ulations  of  the  reflection  of  point  explosion  blast  waves  from  an  inviscid  surface. 
Calculations  were  performed  with  the  REFLECT  code  (Smiley  et  al.,  1981),  and 
the  results  were  intended  to  serve  as  converged  solutions  for  the  ideal-surface  case. 
Figures  25  and  26  present  the  comparisons.  At  low  pressures,  the  REFLECT  code 
results  lie  between  our  smooth-surface  and  ideal-surface  curves  in  the  Mach  reflec¬ 
tion  region  —  especially  at  the  surface  (i.e.,  the  surface  burst  values).  This  indicates 
that  we  have  properly  scaled  the  NP  data.  The  major  difference  appears  in  the  tran¬ 
sition  region.  The  height-of-burst  curves  from  the  REFLECT  code  (denoted  by  the 
dashed  lines)  follow  the  two-shock  theory  all  the  way  to  the  detachment  condi¬ 
tion,  and  then  bulge  upward  and  outward  relative  to  the  NP  curves.  This  trend  is 
prominent  for  reflected  pressures  of  Apn  <  30  psi. 

Figure  27  compares  our  scaled  NP  curves  with  data  from  field  tests  of  large  point 
explosions;  the  typical  yields  were  about  20  KT  (see  Table  2  for  specific  events). 
One  can  see  that  the  data  generally  agree  with  our  curves  in  the  low-pressure  Mach 
reflection  regime  (e.g.,  the  4  psi  curve),  but  the  data  scatter  is  quite  large. 

Figure  28  presents  similar  comparisons  for  intermediate  pressures.  In  a  review  of 
a  preliminary  version  of  our  NP  data,  F.  Sauer  pointed  out  that  the  Tumbler- 
Snapper  series  had  some  high-quality  pressure  measurements  at  the  right  ground 
ranges  and  heights  of  burst  (Event  1  at  HOB  =  754  ft/KT1/3,  Event  2  at  HOB  = 
1009  ft/KT1/3,  and  Event  3  at  HOB  =  1055  ft/KT1/3)  to  clarify  the  shape  of  the 
knees  in  the  transition  region  for  A pr  —  6  to  10  psi.  Subsequently,  H.J.  Carpenter 
has  confirmed  that  the  NOL  measurements  on  these  events  were  indeed  of  high 
quality,  and  he  re- analyzed  the  peak-pressure  versus  range  curves.  Using  these 
data,  we  have  constructed  the  curve  labeled  NOL  in  Figure  28.  At  A p  =  6  psi, 
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Table  2.  Point  exploeion  data  sources. 


SHOB (ft) 

EVENT 

0 

Koa  (HARDTACK  I)1 

0 

Mike  (IVY)1 

2.3 

Butternut  (HARDTACK  I)1 

2.7 

Koon  (CASTLE)1 

3.1 

Sugar  (JANGLE)1 

95 

Harry  (UPSHOT-KNOTHOLE)1 

113 

Annie  (UPSHOT-KNOTHOLE)1 

138 

Met  (TEAPOT)1 

143 

Humboldt  (HARDTACK  II)1 

156 

Apple  II  (TEAPOT)1 

182 

King  (IVY)  1 

196 

Apple  I  (TEAPOT)1 

204 

Priscilla  (PLUMBBOB)1 

205 

Grable  (UPSHOT-KNOTHOLE)1 

205 

Yuma  (REDWING)1 

212 

Galileo  (PLUMBBOB)1 

220 

Wilson  (PLUMBBOB)1 

239 

Morgan  (PLUMBBOB)1 

241 

Bee  (TEAPOT)1 

322 

Climax  (UPSHOT-KNOTHOLE)1 

340 

Hood  (PLUMBBOB)1 

368 

Dog  (TUMBLER-SNAPPER)1 

397 

Easy  (BUSTER)1 

446 

Charlie  (BUSTER)1 

488 

Ruth  (UPSHOT-KNOTHOLE)1 

489 

Dog  (BUSTER)1 

502 

Hamilton  (HARDTACK  II)1 

554 

Franklin  (PLUMBBOB)1 

633 

Charleston  (PLUMBBOB)1 

692 

Wasp  (TEAPOT)1 

702 

Baker  (BUSTER)1 

754 

Able  (TUMBLER-SNAPPER)12 

786 

Encore  (UPSHOT-KNOTHOLE)1 

831 

Rushmore  (HARDTACK  II)1 

1009 

Baker  (TUMBLER-SNAPPER)1'2 

1055 

Charlie  (TUMBLER-SNAPPER)1'2 

1092 

Eddy  (HARDTACK  II)1 

1133 

Mora  (HARDTACK  II)1 

1249 

Lea  (HARDTACK  II)1 

2565 

Dixie  (UPSHOT-KNOTHOLE)1 

Superscript:  1  denotes  Carpenter  (1078)  data  set 


2  denotes  Sauer  (1990)  data  set 
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this  curve  falls  halfway  between  the  NP  curve  and  the  REFLECT  code  curve. 
We  believe  that  this  indicates  that  the  NP  measurements  missed  the  peak  due 
to  gauge-size  effects,  and  the  REFLECT  code  results  overshot  the  peak  due  to  grid 
constraints. 

To  verify  these  statements,  we  have  examined  a  REFLECT  code  calculation  of  a 
HOB  =  600  ft/KT1/3.  The  calculated  pressure-range  curve  is  shown  in  Figure  29. 
The  peak  pressures  decay  with  range  and  agree  with  the  two-shock  theory  for  the 
entire  regular  reflection  region  —  until  the  detachment  angle  ad  is  reached.  Then 
the  pressures  climb  rapidly  until  a  triple  point  (denoted  by  TP)  is  first  observed  in 
the  calculation;  subsequently  the  peaks  decay  rapidly  as  the  Mach  shock  structure 
becomes  better  resolved  on  the  mesh.  This  local  peak  can  dramatically  change 
the  calculated  range-to-effect;  for  example,  the  ground  range  to  A pr  =  30  psi  is 
increased  from  460  ft/KT1/3  to  680  ft/KT1/3  because  of  this  effect.  This  tendency 
for  the  pressure  to  increase  at  transition  is  not  observed  in  any  of  the  data  nor  in 
other  hydrocode  simulations.  Hence,  we  suspect  that  this  local  peak  may  be  an 
artifice  of  the  formulation  of  the  REFLECT  code  calculation. 

The  REFLECT  code  is  a  first-order  Godunov  scheme  that  uses  front  tracking  tech¬ 
niques  for  all  shocks.  The  incident  blast  wave  is  solved  on  a  one-dimensional  grid. 
The  flowfield  inside  the  reflected  shock  is  solved  on  a  two-dimensional  grid  that  con¬ 
forms  to  the  shape  of  the  reflected  shock.  The  two  domains  interact  via  the  jump 
conditions  across  the  reflected  shock.  The  shock  jump  relations  are  also  used  to  ad¬ 
vance  each  local  piece  of  the  reflected  shock.  At  the  reflection  point  a  triangular  cell 
is  used.  The  state  inside  that  cell  is  determined  by  applying  the  two-shock  theory 
(given  the  incident  shock  strength  and  angle).  Hence,  it  is  clear  that  the  method  of 
calculation  forces  the  REFLECT  code  results  to  agree  with  the  two-shock  theory 
throughout  the  regular  reflection  domain.  Near  transition,  the  topology  of  the  grid 
is  changed  to  that  of  a  Mach  reflection  structure.  The  REFLECT  code  specifies 
that  the  state  in  the  quadrilateral  cell  located  at  the  triple  point  is  determined  by 
applying  the  three-shock  theory.  When  the  Mach  stem  is  less  than  one  cell  high 
[i.e.,  for  ground  ranges  from  GR(oj)  to  GR(TP)],  this  method  forces  the  pressure 
to  be  too  large.  Therefore,  we  believe  grid  constraints  in  the  REFLECT  code  could 
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be  the  cause  of  dramatic  increases  in  peak  pressure  at  transition  that  are  not  found 
in  the  data  or  other  (unconstrained)  hydrocode  calculations. 

As  another  illustration,  we  present  the  pressure- vs-range  data  for  Tumbler- 
Snapper  3  (HOB  =  1056  ft/KT1/3)  in  Figure  30.  The  curve  labeled  FS  is 
Fred  Sauer’s  best  estimate  for  this  event  as  predicted  by  his  construction  factor 
method.  The  curve  labeled  NOL  corresponds  to  H.J.  Carpenter’s  evaluation  of  the 
NOL  pressure  records.  The  curve  labeled  NP  comes  from  our  scaled  NP  height- 
of-burst  curves.  It  traces  a  lower  bound  of  all  the  data  from  Tumbler  (where  the 
lower  values  of  peak  pressure  come  typically  from  gauges  with  poorer  response). 
Also  shown  in  Figure  30  are  results  of  an  AMR  calculation  of  a  blast  wave  re¬ 
flection  from  an  inviscid  surface  (Ferguson,  1990).  The  coarse-zoned  calculation 
(Amin  =  10  ft/KT1/3)  gave  peak  pressures  in  agreement  with  the  NP  curve  (gauge 
size  ~  8  ft/KT1/3),  while  the  fine-zoned  calculation  (Amj„  =  1  ft/KT1/3)  gave  peak 
pressures  that  were  consistent  with  the  NOL  curve  (gauge  size  ~  0.3  ft/KT1/3). 
We  infer  from  this  that  the  gauges  on  the  NP  tests  were  too  large  in  diameter  to 
accurately  measure  these  peaks  at  transition. 

Figure  31  depicts  the  entropy  contours  at  various  times  from  the  aforementioned 
fine-zoned  AMR  calculation  of  the  reflection  of  a  point-explosion  blast  wave  from 
am  inviscid  surface  (HOB  =  1056  ft/KT1/3  =  322  m/KT1/3).  Entropy  was  selected 
as  the  display  variable  because  it  is  the  most  sensitive  indicator  of  a  contact  surface 
emanating  from  a  triple  point.  The  contour  plots  show  the  transition  from  a  regular 
reflection  shock  structure  to  a  Mach  reflection  shock  structure.  An  entropy  layer 
begins  to  appear  along  the  wall  at  a  GR  ~  280  m  =  919  ft  ( a  =  41°)  and  is  clearly 
defined  at  a  GR  =  388  m  =  1273  ft  (a  =  50°).  Eventually  the  entropy  layer  lifts 
off  the  wall  and  an  identifiable  Mach  stem  is  observed,  e.g.,  GR  cs  630  m  =  2067 
ft  (a  =  63°).  If  one  defines  transition  as  the  formation  of  an  entropy  layer  near 
the  surface,  then  this  calculation  indicates  that  transition  occurs  at  an  a  =  41°  to 
50°,  while  the  theoretical  detachment  angle  for  this  case  is  a<j  =  57°.  Hence,  this 
inviscid  AMR  calculation  indicates  that  the  shock  structure  transitions  well  before 
the  detachment  condition  is  reached. 
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As  the  Mach  stem  evolves,  the  gradient  in  entropy  from  the  triple-point  region  (i.e., 
the  contact  surface)  weakens  and  then  disappears,  e.g.,  GR  =  1052  m  =  3451  ft 
(a  =  73°).  The  shock  structure  appears  to  be  similar  to  the  Neumann  Reflection 
NR  structure  described  by  Henderson  and  Colella  (1990)  for  weak  shock  reflections 
from  wedges.  Recall  that  for  a  NR  structure,  the  outer  shock  is  continuously  curved 
(i.e.,  there  is  no  kink  at  the  triple  point)  and  the  reflected  shock  spreads  out  into 
a  dispersed  wave  near  the  front.  Thus,  the  triple  point  becomes  a  dispersed  wave 
region  devoid  of  any  contact  surface  or  slipline.  Notice  also  in  Figure  31,  that  the 
gradients  behind  the  reflection  point  are  much  steeper  than  the  gradients  in  the 
incident  blast  wave.  Thus,  one  needs  a  much  finer  mesh  or  much  smaller  gauges  to 
accurately  capture  peaks  for  such  Neumann  Reflections  (compared  to  what  would 
be  needed  to  resolve  the  incident  blast  wave). 

To  summarize,  the  above  comparisons  suggest  that  the  REFLECT  code  results 
overshoot  at  transition  due  to  grid  constraints.  The  NP  experiments  missed  the 
peak  at  transition  because  of  gauge-size  effects;  hence,  smaller  gauges  (<  1  ft/KT1/3 
=  0.36  mm/. 5  g1/3)  are  needed  to  accurately  measure  such  peaks.  Fortunately,  PCB 
has  recently  developed  smaller  pressure  gauges  that  may  significantly  improve  the 
peak  pressure  measurements  (see  Section  9).  The  AMR  results  were  consistent  with 
the  NOL  data  and  show  that  although  the  overpressure  does  not  locally  peak  near 
transition,  the  pressure-range  curve  does  flatten  and  lead  to  more  pronounced  knees 
in  the  low  pressure  HOB  curves  than  indicated  by  the  NP  data.  These  calculations 
also  showed  that  the  shock  structure  transitions  to  a  dispersed  wave  system  in  the 
low-pressure  regime. 
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SECTION  8 
CONCLUSIONS 


1.  NP  Height-of- Burst  Curves. 

Measurements  made  on  the  0.5-g  NP  experiments  form  a  very  extensive  and 
self-consistent  height-of-burst  data  set.  Because  the  NP  data  is  self-consistent 
(i.e.,  the  same  charge,  the  same  gauges  and  the  same  surface  were  used),  one 
can  detect  subtle  effects  such  as  boundary  layer  effects. 

The  shape  of  the  NP  height-of-burst  curves  agrees  with  data  from  height-of- 
burst  experiments  performed  with  much  larger  HE  charges.  The  NP  data  are 
also  consistent  with  the  low  fidelity  (i.e.,  poor  response)  point  explosion  data. 
Nevertheless,  comparison  with  high  fidelity  point  explosion  data  (e.g.,  the  NOL 
data)  and  recent  AMR  hydrocode  calculations  has  demonstrated  that  the  NP 
peak  pressures  Eire  about  20%  too  low  in  the  low-pressure  transition  region  — 
due  to  gauge-size  limitations.  Hence,  smaller  gauges  (or  larger  explosions)  are 
needed  to  more-accurately  define  the  peak  pressures  in  the  transition  region. 

2.  NR  Shock  Structure. 


According  to  recent  hydrocode  simulations,  the  reflected  shock  structure  transi¬ 
tions  to  a  Neumann  Reflection  (NR) — not  a  classical  three-shock  Mach 
structure — for  weak  shock  reflections.  This  is  confirmed  by  photography  (in 
Fig.  5).  There  is  a  strong  rarefaction  wave  associated  with  this  NR  structure 
which  causes  a  steep  pressure  gradient  behind  the  foot  of  the  Mach  stem  and 
makes  it  difficult  to  resolve  the  peak  (experimentally  or  computationally). 

3.  Transition  Criterion. 

The  transition  criterion  is  not  known  for  weak  blast  waves.  According  to  pressure 
measurements  on  HE  experiments  (Figs.  9,  10  and  11)  and  point  explosion  tests, 
the  transition  seems  to  occur  (i.e.,  reflected  pressures  start  to  deviate  from  the 
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two-shock  theory)  before  the  detachment  angle  is  reached:  oct  —  40° (A pn  =  30 
psi)  and  ar  —  28°(  ApK  =  8  psi).  It  is  r.  jt  known  whether  this  is  a  gasdynamic  effect 
or  a  boundary  layer  effect,  or  whether  it  is  caused  by  measurement  inaccuracies. 

4.  Boundary  Layer  Effects. 

Boundary  layer  effects  on  peak  reflected  pressures  were  estimated  from  the  data 
taken  on  the  hydrodynamically-smooth  surface.  Negligible  effects  were  found  at 
small  ground  ranges.  However,  boundary  layer  effects  probably  caused  a  deficit 
in  range  to  the  same  overpressure  that  increased  logarithmically  with  ground 
range: 

6(%)  =  7.2  +  14.6  log(Gf2) 

for  peak  pressures  below  15  psi.  The  effects  of  the  boundary  layer  on  the  tran¬ 
sition  ground  range  could  not  be  evaluated;  comparisons  with  inviscid-surface 
hydrocode  calculations  are  required  for  that  purpose. 

5.  Other  Real- Surface  Effects. 

Rough- wall  boundary  layer  effects  were  evaluated  from  the  pressure  data  taken 
on  the  #40  sandpaper  surface  by  comparison  with  the  smooth-surface  results. 
No  effect  on  the  peak  pressures  was  found  in  the  regular  reflection  region.  A 
slight  decrease  in  the  ground  range  for  transition  was  observed.  After  transition, 
a  small  but  systematic  depression  of  the  peak  pressures  was  measured  in  the  knee 
region  of  the  curves.  Below  the  knee  region,  the  large  roughness  had  essentially 
no  effect  on  the  peak  pressure  HOB  curves. 

The  effects  of  mass  loss  due  to  surface  porosity  were  evaluated  from  the  pressure 
data  taken  on  the  Filtrokelit  surface.  Porosity  caused  no  effect  on  the  peak 
pressures  in  the  regular  reflection  region.  Porosity  had  some  effect  on  the  ground 
range  to  transition,  and  this  effect  increased  at  lower  pressures  and  higher  burst 
heights.  The  large  porosity  employed  in  these  tests  dramatically  modified  the 
height-of-burst  curves  from  transition  to  the  knee  region,  but  below  the  knee 
region,  porosity  had  no  effect  on  the  curves. 
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The  effects  of  surface  compressibility  were  evaluated  from  the  pressure  data 
taken  over  a  deep  snow  layer  (Wisotski,  1990).  Dramatic  changes  to  all  portions 
of  the  height-of-burst  curves  were  observed  for  all  reported  pressures  (A pn  =  3 
to  30  psi).  This  may  be  a  deep-snow  effect  that  would  be  reduced  if  a  scaled 
snow  depth  is  used. 
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SECTION  9 

RECOMMENDATIONS 


We  offer  the  following  recommendations  for  improvements  to  the  experimental  tech¬ 
nique,  additional  tests  to  explore  real-surface  effects,  and  supporting  hydrocode 
calculations. 

1.  Smaller  Pressure  Gauges. 

Smaller-diameter  pressure  gauges*  all  needed  to  accurately  measure  the  very 
peaked  waveforms  in  the  low-pressure  transition  region  and  to  investigate  the 
double-peaked  waveforms  in  the  double-Mach  reflection  in  the  high-pressure 
regime.  First,  shock  tube  tests  should  be  performed  with  such  st^all  gauges 
mounted  next  to  the  Kistler  603B  gauges.  This  will  allow  direct  comparisons 
of  the  recorded  pressures  measured  at  the  same  position  on  the  same  test,  by 
which  we  can  evaluate  potential  improvements  in  measured  rise  times  and  peak 
values.  Side-on  tests,  end -wall  tests  and  wedge-reflection  tests  are  also  needed. 
Then,  the  smaller  gauges  can  be  used  in  future  HOB  tests. 

2.  Larger  Explosion  Chamber. 

Currently  the  HOB  curves  for  A p  =  1,  1.5  and  2  psi  are  based  on  extrapolations 
of  the  pressure- range  curves.  A  larger  explosion  chamber  (~  2m  by  2m  by  1.5m) 
should  be  built  so  that  the  low-pressure  HOB  curves  can  be  constructed  from 
measurements — not  extrapolations. 

3.  Photography  of  Mach  Reflection  Transition. 

Shadow- schlieren  photography  should  be  used  to  systematically  investigate  the 
shock  structure  in  the  transition  region.  This  can  be  used  to  estimate  Mach 
transition. 


*For  example,  the  PCB  Corporation  has  recently  developed  crystal  gauges  with 
a  sensitive  element  of  0.015  inches  across,  which  have  an  output  of  hundreds  of 
millivolts  per  psi. 
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4.  High-Pressure  HOB  Curves. 

Additional  tests  with  intermediate  HOBs  and  ground  ranges  are  needed  to  more- 
fully  define  the  NP  height-of-burst  curves  in  the  high-pressure  region  (A p  =  100 
to  1000  psi). 

5.  Real-Surface  Effects. 

The  height-of-burst  tests  should  be  continued  to  produce  a  complete  set  of  NP 
height-of-burst  curves  (i.e.,  down  to  pressures  of  A pr  =  1  psi)  for  both  hydro- 
dynamically  smooth  and  rough  surfaces. 

In  addition,  height-of-burst  tests  should  be  performed  over  other  nonideal  sur¬ 
faces,  for  example:  Astro-Turf  carpet  could  be  used  to  simulate  shock  diffraction 
effects  caused  by  forests  and  vegetation,  polyurethane  foam  could  be  used  to  sim¬ 
ulate  the  compressibility  effects  of  snow-covered  surfaces,  and  Freon  layers  could 
be  used  to  simulate  density  effects  in  boundary  layers. 

6.  Hydrocode  Simulations. 

Inviscid  hydrocode  simulations  of  selected  height-of-burst  tests  would  be  use¬ 
ful  to  examine  the  transition  process  for  weak  blast  waves,  and  to  completely 
eliminate  wall  boundary  layer  effects  (e.g.,  by  using  an  inviscid  slip  boundary 
condition  on  the  surface).  A  second-order  code  with  adaptive  mesh  refinement 
and  front  tracking  is  required  for  such  studies. 

7.  Field  Tests. 

As  already  mentioned,  the  pressure  waveforms  are  extremely  peaked  in  the  low- 
pressure  transition  region;  this  makes  it  difficult  to  accurately  measure  the  peak 
pressure  values  in  this  region.  In  item  1  above,  we  recommended  using  the 
smaller-size  PCB  gauges  to  increase  the  temporal  resolution  of  the  pressure 
measurements  in  the  0.5-g  charge  height-  of-burst  tests.  To  increase  the  reso¬ 
lution  even  further,  one  must  increase  the  scale  of  the  explosion.  For  example, 
0.5- kg  charges  (which  could  be  used  at  the  EMI  test  site  at  Winters weiler)  would 
increase  the  explosion  scale  by  a  factor  of  10,  and  8-lb  charges  (which  could  be 
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used  at  the  DNA  test  site)  would  increase  the  explosion  scale  by  a  factor  of  19.  It 
would  be  useful  to  perform  a  few  field  tests  with  such  larger  charges  and  small 
PCB  gauges  —  to  check  the  waveforms  and  peak  pressures  in  the  transition 
region  for  burst  heights  around  1000  ft/KT1/3. 
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Figure  2.  Construction  of  NP  charges:  (a)  spherical  charge;  (b)  hemispherical  charge; 
(c)  mounting  bracket  for  the  hemishperical  charge. 
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Figure  3.  Schematic  of  the  photography  system. 
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Figure  4.  Schematic  of  the  pressure  gauge  mount. 
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Figure  5. 


Shadow  photographs  showing  the  blast  reflection  over  the  Makrolon 
surface  for  HOB  =  303  mm  =  3.8  m/kg1/3:  (a)  regular  reflection 
structure  near  transition  ( GR  =  4./3  m/kg1/3,  o  =  51.2  , 

Api  =4.1  psi  and  Apr  cs  10  psi);  (b)  Neumann  reflection 
shock  structure  well  past  transition  (GR  =  9.02  m/kg1/3, 
a  =  67°,  Api  =  2  psi  and  Apr  ~  5  psi). 
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Figure  6.  Pressure- range  curves  for  0.5-g  NP  charges  (subscripts  I 
and  R  denote  incident  and  reflected  values,  respectively). 


pressure  height-of-burst  curves  for  a  hydrodynamically- smooth  surface 
e  intermediate-pressure  regime:  notation  is  the  same  as  Figure  7. 


Figure  9.  Peak-pressure  height-of-burst  curves  for  a  hydrodynamically- smooth  surface 
in  the  high-pressure  regime:  notation  is  the  same  as  Figure  7. 


Carpenter,  1978 


Comparison  of  the  NP  height-of  burst  curves  for  a  smooth  surface 
with  data  from  450-kg  TNT  charges  (Ap  =1  to  15  psi). 
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Figure  12.  Comparison  of  the  NP  height-of- burst  curves  for  a  smooth  surface 
with  data  from  1-kg  RDX  charges  (symbol  x  ). 
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Figure  13.  Comparison  of  the  NP  data  with  the  height-of-burst  curves  for 
3.6-kg  PBX-9404  charges  (Carpenter  and  Brode,  1974). 
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Figure  11.  Peak-pressure  height-of-burst  curves  for  a  smooth  surface  (symbol  0),  a  rough 
surface  (symbol  R)  and  a  porous  surface  (symbol  P). 
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Figure  IS.  Peak-pressure  decay  with  range  for  surface  bursts: 
circles  denote  an  ideal  surface;  squares  denote  an 
extrapolated  surface  burst 
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Figure  16.  Ground  range  decrement,  6,  versus  ideal  range. 
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Figure  17.  Peak-pressure  height-of-burst  curves  for  a  smooth  surface  (symbol  0) 
and  an  ideal  surface  (symbol  /)  in  the  intermediate-pressure  regime. 
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Figure  18.  Peak-pressure  height-of-burst  curves  for  a  smooth  surface  (symbol  0)  and  an  ideal  surface  (symbol  /) 
in  the  low-pressure  regime.  The  curves  for  A p  =  2,  1.5  and  1  psi  are  extrapolated. 


arision  of  reflected  pressures  on  a  deep  snow  surface  (Wisotski, 
with  the  smooth -surface  NP  results  (low-pressure  regirr\e). 


curves  for  a  smooth  surface  (symbol  O)  and  an  ideal  surface 
o  a  1-KT  point  explosion  by  function  SF  for  Ap  =  1  to  15  psi 


REFLECT  code  calculations 


GR  ( ft/KT1/J) 

Comparison  of  the  scaled  NP  height-of-burst  curves  with  REFLECT  code  calculations 
of  point  explosions  (Smiley  et  al.,  1981)  for  Ap  =  1  to  15  psi. 


REFLECT  code  calculations 


GR  ( ft/KT  J) 

Comparison  of  the  scaled  NP  height-of-burst  curves  with  REFLECT  code  calculations 
of  point  explosions  (Smiley  et  al.,  1981)  for  Ap  =  6  to  200  psi. 


Point  explosion  data 


Comparison  of  the  scaled  NP  height-of-burst  curves  with  point  explosion  data, 


Carpenter  (1978) 


GR  {fl/KT1/J) 

Comparison  of  the  scaled  NP  height-of-burst  curves  with  point  explosion  data, 


A  PR  (psi) 


200  400  600  800  1000 


GH(ft/KT1/3) 


Figure  29.  Peak  pressure  vs.  range  for  a  point  explosion  at  HOB  =  600  ft/KTi# 
from  the  REFLECT  code  calculations  (Smiley  et  al.,  1981). 
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Figure  30.  Comparison  of  peak-overpressure  vs.  range  curves  for  HOB  =  1056  ft/KT1^ 


Figure  3 1 .  Entropy  contours  showing  the  reflection  of  a  point-source  blast  wave 
from  a  plane;  HOB  =  1056/KT1/3  *  322  m/KT1'3  (Continued). 


HOB  (m/kgi/3) 


HOB  (m/kgi/3) 


GR  (m/kgi/3) 


Figure  A3.  Comparison  of  the  NP  height-of-burst  curves  for  shock  arrival  time,  TOA, 
for  a  smooth  surface  vs.  a  porous  surface. 
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Figure  A4. 


GR  (m/kgi/3) 

Comparison  of  the  NP  height-of-burst  curves  for  positive  phase  duration,  x. 
for  a  smooth  surface  vs.  a  porous  surface. 
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